Derealization of two-particle ring near the Fermi level of 2d Anderson model 
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We study analytically and numerically the problem of two particles with a long range attractive 
interaction on a two-dimensional (2d) lattice with disorder. It is shown that below some critical 
disorder the interaction creates delocalized coupled states near the Fermi level. These states appear 
inside well localized noninteracting phase and have a form of two-particle ring which diffusively 
propagates over the lattice. 

PACS numbers: 72.15.Rn, 71.30+h, 74.20.-z, 05.45.Mt 



Recently a great deal of attention has been attracted to 
the problem of interaction effects in disordered systems 
with Anderson localization 10,0. From the theoretical 
point of view the problem is rather nontrivial. Indeed, 
even if a great progress has been reached in the theo- 
retical investigation of the properties of localized eigen- 
states H still the analytical expressions for interaction 
matrix elements between localized states arc lacking. In 
spite of these theoretical difficulties it has been shown re- 
cently that a repulsive or attractive interaction between 
particles can destroy localization and lead to a propaga- 
tion of pairs in the noninteracting localized phase. This 
two interacting particles (TIP) effect has been studied 
recently by different groups 0-11 and it has been un- 
derstood that the derealization of TIP pairs is related 
to the enhancement of interaction in systems with com- 
plex, chaotic eigenstates. Such an enhancement had been 
already known for parity violation induced by the weak 
interaction in heavy nuclei jl2| where the interaction is 
typically increased by a factor of thousand. However, 
since there the two-body interaction is really weak the 
final result still remains small. On the contrary, for TIP 
pairs in the localized phase the enhancement of inter- 
action qualitatively changes the dynamics leading to a 
coherent propagation of TIP on a distance l c being much 
larger than the pair size and one-particle localization 
length l\ . The enhancement factor n is determined by the 
density of two-particle states p2, coupled by interaction, 
and the interaction induced transition rate T2 between 
noninteracting eigenstates, so that k = T2P2- At K ~ 1 
the interaction matrix element becomes comparable with 
two-particle level spacing and the Anderson localization 
starts to be destroyed by interaction. For excited states 
the TIP density P2 is significantly larger than the one- 
particle density p and the derealization can be reached 
for relatively weak interaction if l\ is large. However, 
when the excitation energy e above the Fermi level de- 
creases then p2 becomes smaller and it approaches the 
one-particle density p at low energy: p2 ~ ep 2 ||. As a 
result the value of k also drops with e so that the dereal- 
ization of TIP pairs practically disappears near the Fermi 
energy. This result has been found in |5|,|l3[] in the ap- 



proximation of the frozen Fermi sea created by fermions. 
Recent numerical studies of TIP pairs with short range 
interaction near the Fermi level [fl4| confirmed these the- 
oretical expectations. 
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FIG. 1. Probability distributions / and f d for TIP in 2d 
disordered lattice of size L — 40, and interaction of radius 
R — 12 and width AR — 1. Left column, one-particle prob- 
ability / for W = 8V: ground state at U — (top); ground 
state with binding energy AE = —1.05V at U = —2V (mid- 
dle); coupled state with AE = -0.19V at U = -2V (bottom). 
Right column: / for coupled state, compare to bottom left, 
at W = 12V and U = -2V with AE » -0.191/ (top); in- 
ter-particle distance probability fd related to the middle left 
case (middle); fd related to the bottom left case (bottom). 
All data are shown for the same disorder realisation. Color 
2d density plots for these data are given in Fig. Ibis of Ap- 
pendix. 
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In this paper we discuss another type of situation in 
which TIP derealization (see Fig. 1) takes place mainly 
due to geometrical reasons and not due to the relation 
p 2 ^> p. As a result the TIP pair can be delocalized in 
a close vicinity to the Fermi level that opens new inter- 
esting possibilities for interaction induced derealization 
in the localized noninteracting phase. To study this new 
situation we investigate a model with a long range at- 
tractive interaction between particles in the 2d Anderson 
model. In this case the particles can rotate around their 
center of mass, being far from each other and keeping 
the same energy, while the center can move randomly in 
two dimensions. As a result the system has effectively 
three degrees of freedom that makes it rather similar to 
the case of one particle in the 3d Anderson model where 
derealization takes place at sufficiently weak disorder. 
A somewhat similar situation has been studied recently 
for particles with Coulomb interaction but only excited 
states were considered there and the derealization was 
attributed to the large ratio pi I 'p @|. Here we show 
that in fact the conditions for derealization are much 
less restrictive. 

To illustrate the above ideas let us first discuss the case 
of only two particles with attractive interaction U(r) < 
in the 2d Anderson model described by the Schrodinger 
equation 



(E ni +E n2 +U (m - n 2 ))Vm, 



-l,n 2 M\ 
^ni ,n 2 +l + ^111,112 — 1 ni,n 2 ■ 

Here ni2 are the indices of the two particles on the 2d 
lattice with L 2 sites and periodic boundary conditions, V 
is the hopping between nearby sites and the random on- 
site one-particle energies E ni 2 are homogeneously dis- 
tributed in the interval [—W/2,W/2]. The long range 
attractive interaction depends only on the distance be- 
tween particles r = ||ni — nail and is equal to a constant 
U < if |r — R\ < AR and zero otherwise. The value 
of r is determined as the minimal inter-particle distance 
on the periodic lattice. Thus the interaction takes place 
only inside a ring of radius R and width AR, and we as- 
sume that R ^> AR > 1. For U = the eigenstates are 
given by the product of two one-particle (noninteracting) 
eigenstates which are always localized in 2d in a presence 
of disorder |16| . 

In the limit ofvery strong attractive interaction |f7| 3> 
V the TIP coupled states form the energy band of width 
~ 16V around E ~ —\U\ (we consider only the states 
symmetric in respect to particle interchange). For states 
in this band the particles are located always inside the 
ring which center can move over the 2d lattice. Since 
\U\ ^> V these states are decoupled from all other 
states with particles outside the ring. In the ring the 
Schrodinger equation is in fact rather similar to the case 
of 3d Anderson model of one particle. In this analogy 
the number of sites inside the ring Mr « 2irRAR deter- 
mines the effective number of 2d planes placed one over 



another in the third z-dimension (length size L z = Mr). 
In this 3d model the effective strength of disorder is ap- 
proximately 2W since the diagonal term is now the sum 
of two E n values. Also one site is coupled with Z = 8 
neighbours contrary to Z = 6 for 3d case (assuming 
AR ^> I). Since in 3d the Anderson transition at the 
band center takes place at W c = 2.75ZV = 16.5V [0, 
we expect that TIP states inside the ring will be delo- 
calized in the middle of the band when 2W/ZV = 2.75 
that gives the transition at W c2 ~ 11V. This estimate 
is in agreement with numerical simulations of the model 
(0) jlSj ]- Of course, since the size in the third direction 
is finite the eigenstates will be eventually localized. But 
their localization length l c will make a sharp jump from 
l c - I at W > W c2 to l c ~ exp(g) > I at W < W c2 that 
follows from the standard scaling theory in 2d | Ml|[l5j . 
Here g is the conductance of the quasi-two-dimensional 
layer of width L z = Mr. As usual g = E c /A\ where 
E c = D/L 2 is the Thouless energy, Ai ~ V/(L 2 L Z ) is 
the level spacing H and the diffusion rate in the lattice 
model is D ~ V(V/W) 2 . As a result for W < W c2 the 
TIP derealization length jumps to exponentially large 
value l c ~ exp(2irRAR(W c2 /W) 2 ). In these estimates 
we assumed that l\ > AR > 1 since if Ai? l\ the ma- 
jority of states inside the ring are noninteracting and can 
be presented as the product of one-particle eigenstates. 
We also note that for W < W c2 there is an energy interval 
around the band center with delocalized states where the 
TIP ring diffuses with the rate D 2 ~ V{W c2 /V) 2 . When 
W decreases the mobility edge approaches the bottom of 
the band as it happens in 3d Anderson model. 

The above arguments presented for the case \U\ ^> V 
indicate that it is possible to have a similar TIP dereal- 
ization at moderate value of U ~ V near the Fermi level. 
To investigate this case we rewrite the equation (^) in 
the basis of the noninteracting eigenstates that gives 



ni2 .m' 2 ^(■ni. ,m' 2 



(2) 



Here Xm 1 ,m 2 are eigenfunctions of the TIP problem 
written in the basis of one-particle eigenstates <p m 
with eigenenergies Em. The matrix UQ > < 

represents the two-body matrix elements of interac- 
tion U{\\\ — n 2 ) between noninteracting eigenstates 
\4>m 1 4>m 2 ) an d \<f> m ' 4> m > )■ The Fermi sea is determined by 

the restriction of the summation in @ to m\ \ > with 
energies E ( / } > Ep, where Ep is the Fermi energy re- 

m l,2 

lated to the filling factor p. We choose the case with half 
filling p — 1/2 for which Ep ss 0. In this way our model 
corresponds to the approximation of frozen Fermi sea suc- 
cessfully used for the Cooper problem 19 1 . As it was done 
by Cooper we also introduce the high energy cut-off de- 
fined by the condition 1 < m 1 + m 2 < M. This rule 
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determines an effective phonon frequency uj d oc M/L 2 . 
We fix a = L 2 /M « 15 since uj d should be independent 
of the system size L [|o) . We checked that the results are 
not affected by a variation of a in few times. The first 
studies of the TIP model with frozen Fermi sea was done 
by Imry || with the aim to take into account the effect 
of finite fermionic density and then was also analyzed in 
p3[ . Recently a similar model was investigated for the 
case of Hubbard attraction in 3d jl4j . 
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FIG. 2. Dependence of IPR £ on the binding energy AE /V 
for U = -2V, W = 8V" and AR = 1: L — 20 (o), L — 22 
(triangle down), L = 26 (diamond), L = 30 (triangle up), 
L = 40 (square); full/empty symbols are for R = 8/R — 12; 
the shaded band shows the variation of £ at £f and U = for 
20 < £ < 40. Insert shows dependence of £ on disorder VK/ V 
at R = 8 for the ground state at U = — 21/, L = 30 (x) and 
at U = 0, L = 30 (+), and for the states at the derealization 
border with binding energy AE ~ AE C at U = —2V (*). 

To study the eigenstate properties of our model we di- 
agonalize numerically the Hamiltonian (|^) and rewrite 
the eigenfunctions in the original lattice basis. In this 
way we determine the two-particle probability distribu- 
tion F(ni,ri2) from which we extract the one particle 
probability /(ni) = X)n 2 ^( ni ' n2 ) an< ^ ^he probabil- 
ity of inter-particle distance /d(r) = ^ n2 + n2 ' n2 ) 
with r = ni ri2. The binding energy of an eigenstate 
in (§) is AE = E - 2E F w E since E F » 0. For 
the ground state with energy E g the coupling energy is 
A = 2Ep — E g . The typical examples of probability dis- 
tributions are shown in Fig. 1. They clearly show that 
the ground state in the presence of interaction remains 
localized and the particles stay on distance R from each 
other. However, there are states with negative binding 
energy (AE < 0) which are delocalized by interaction 
and for which the particles move around the ring in agree- 



ment with discussion of model (|l|) at \U\ 3> V. We stress 
that this derealization of coupled states (AE < 0) takes 
place in the well localized one-particle phase. However, 
at very strong disorder this derealization disappears (see 
top right case in Fig. 1). 

To analyze the derealization of states with negative 
binding energy AE in a more quantitative way we de- 
termine the inverse participating ratio (IPR) £ for one- 
particle probability l/£ = (X)n/ 2 ( n ))> w here brackets 
mark the averaging over 100 disorder realisations. In this 
way £ gives the number of lattice sites occupied by one 
particle in an eigenstate. The dependence of £ on AE 
and W is shown in Fig. 2 for different lattice sizes L in 
the presence of interaction. This figure shows that near 
the ground state the interaction creates states which are 
even more localized than in the absence of interaction (£ 
is significantly smaller than at U = 0, see insert Fig. 2). 
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FIG. 3. Level spacing statistics P(s) for TIP coupled states 
at L = 30, U = -2V, W = 8V, R = 8, AR = 1 in the lo- 
calized phase near the ground state inside the energy interval 
—A < AE < — 3A/4 (o) and in the delocalized phase for en- 
ergies inside -A/4 < AE < (•); here AE C fa -0.3A and 
the statistics is done over 3000 disorder realisations. Full lines 
show the Poisson distribution and the Wigner surmise. 

In fact for -A < AE < AE C < the IPR value even 
slightly drops with the increase of L. However for the 
states with binding energy AE C < AE < the situ- 
ation becomes different and £ grows significantly with 
L while the change of IPR at U = with L is rather 
weak (see shaded band in Fig. 2). The critical value of 
the binding energy AE C can be defined as such an en- 
ergy at which £ remains independent of L. In this way 
AE C determines the mobility edge for coupled states so 
that at given U and W the TIP eigenstates are local- 
ized for -A < AE < AE C while for AE C < AE < 
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the states becomes delocalized (see an example in Fig. 
1). In agreement with this picture £ varies up to 30 
times when AE changes from —A up to 0. This vari- 
ation grows with L and the interaction radius R since 
the system becomes more close to the effective 3d An- 
derson model as it was discussed above. The qualitative 
change of the structure of the eigenstates leads also to a 
change in the level spacing statistics P(s) (Fig. 3). Near 
the ground state the statistics is close to the Poisson dis- 
tribution Pp(s) = exp(— s) typical for the localized An- 
derson phase pl|] while for AE C < AE < it approaches 
to the Wigner surmise Pw( s ) = 7rs exp(— 7rs 2 /4) /2 cor- 
responding to the delocalized phase fl2j]]. 
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FIG. 4. Dependence of the derealization border AE C and 
the binding energy of the ground state —A on disorder W/V. 
The values of AE C are shown by (o)/(») for U = -V / - 2V; 
the values of —A are shown by upper/lower dashed line for 
U = — V I — IV and L = 30 (vertical intervals give variation 
for 20 < L < 30). 

The variation of the derealization border AE C for TIP 
coupled states with disorder strength and interaction is 
shown in Fig. 4. While the coupling energy A grows with 
U and W, the mobility edge AE C < 0, on the contrary, 
disappears at strong W. According to the data of Fig. 4 
all states with binding energy AE < become localized 
for W > W c2 ~ 9.5V (U = -2V) and W > W c2 ~ 8V 
(U = —V). This shows that at weaker interaction a 
weaker disorder is required to have delocalized coupled 
states. As it follows from Fig. 4, at small disorder W 
the derealization border AE C becomes closer and closer 
to the ground state. This means that at weak disor- 
der the derealization will take place for excited states 
with low energy. For W <C W C 2 and U ~ — V the dif- 
fusion rate of delocalized TIP ring can be estimated as 
D 2 ~ V(W C 2/W) 2 p2| . Further studies are required to 
determine the dependence of W C 2 on W at \U\ <C V . 



In conclusion, our results show that long range at- 
tractive interaction between two particles in 2d leads to 
the appearance of delocalized diffusive states near the 
Fermi level inside the well localized noninteracting phase. 
It would be interesting to understand what will be the 
consequences of this derealization for real many-body 
fcrmionic problem with attractive interaction. It is pos- 
sible that obtained results will be also relevant for elec- 
trons with Coulomb repulsion. Indeed, in this case at 
very weak disorder each electron oscillates near an equi- 
librium position and the two-body interaction can be con- 
sidered as an effective harmonic attraction |23| . 
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FIG. Ibis. Color 2d density plots for the data of Fig.l with the same ordering of figures. Blue corresponds to the 
minimum of the probabilty distribution and red to the maximum. The first four figures are drawn in logarithmic scale 
while two figures at the bottom arc in linear scale. Blue/Red color corresponds to: / = 1.3 x 10 -10 // = 0.2 (top 
left), / = 1.1 X 10- g /f = 0.26 (middle left), / = 1.5 X 10~ 6 // = 0.014 (bottom left), / = 1.14 x lO" 11 // = 0.073 
(top right), f d = 3.8 x 10- 7 // d = 0.1 (middle right), f d = 1.5 x 10 _5 //<i = 0.0032 (bottom right). 



